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Abstract 

A theory of fractional kinetics of glial cancer cells is presented. A 
role of the migration-proliferation dichotomy in the fractional cancer cell 
dynamics in the outer-invasive zone is discussed an explained in the frame¬ 
work of a continuous time random walk. The main suggested model 
is based on a construction of a 3D comb model, where the migration- 
proliferation dichotomy becomes naturally apparent and the outer-invasive 
zone of glioma cancer is considered as a fractal composite with a fractal 
dimension Df r < 3. 

KEYWORDS: Glioma; Migration-Proliferation Dichotomy; Fractional 
Kinetics. 


1 Introduction 

Brain tumors result from the uncontrolled growth of abnormal cells, destruction 
of normal tissues, and invasion of vital organs. These processes can be subdi¬ 
vided into many types based on several classification characteristics and involve 
any of the cell types found in the brain, such as neurons, glial cells, astrocytes, 
or cells of the meninges mm- The mechanisms behind cancer progression result 
from the accumulation of one or a few specific mutations that disrupt biological 
pathways like growth factor signaling, DNA damage repair, cell cycle, apoptosis 
and cellular adhesion [3]. Among all possible cancer cell genotypes, leading to 
six main alternations of malignant growth [3] , cell motility and invasion are the 
most important for the present consideration. 

Glioma is one of the most recalcitrant brain disease, with an optimal therapy 
treatment survival period of 15 month and most tumors recur within 9 months 
of initial treatment HE]. One of the main possible reason of such devastat¬ 
ing manifestation is the migration proliferation dichotomy of cancer cells. This 


1 


phenomenon has been firstly observed at clinical investigations 0 m, where it 
has been shown that in the outer invasive zone glioma cancer cells possesses a 
property of high motility, while the proliferation rate of these migratory cells is 
essentially lower than in the tumor core. This anti-correlation between prolifera¬ 
tion and migration of cancer cells, also known as the Go or Grow hypothesis (see 
discussions in 00), suggests that cell division and cell migration are temporally 
exclusive phenotypes 0. The phenomenon that tumor cells defer proliferation 
for cell migration was also experimentally demonstrate^ in m\mm- The 
switching process between these two phenotypes is still not well understood. 
Moreover, it should be mentioned that conflicting data appear in the literature 
concerning the Go or Grow hypothesis; details of discussions on this can be 
found in 00. 

Extensive theoretical modelling follow this finding. A switching process be¬ 
tween these two phenotypes still is not well understood, and many efforts are 
directed to develop relevant models with relevant mechanisms of switching of 
the glioma cells, resulting in several phenomenological models. Comprehensive 
discussions of these models one can find in ® nsa nsi Ha ns- It was sug¬ 
gested by Khain et al ® that the motility of cancer cells is a function 
of their density. Multi-parametric modelling of the phenotype switching was 
considered in [221. The agent-based approach to simulate multi-scale glioma 
growth and invasion was used in EU[22]. Subdiffusive cancer development on a 
comb was studied in [231 . where a continuous time random walk (CTRW) was 
firstly suggested for metastatic cancer development. A stochastic approach for 
the proliferation-migration switching involving only two parameters was pro¬ 
posed in [24) :25j where the transport of cancer cells was formulated in terms of 
the CTRW, as well. ‘Go or Grow’ mechanism was proposed in [15], where the 
transition to invasive tumor phenotypes can be explained on the basis of the 
oxygen shortage in the environment of a growing tumor. Phenotypic switching 
due to density effect was also suggested in [MIlS. Both numerical and an¬ 
alytical approaches were developed in mi to study the glioma propagation in 
the framework of reaction-diffusion equations, where the phenotype switching 
depends on oxygen in a threshold manner. Collective behavior of brain tumor 
cells under the hypoxia condition was studied in m- 

A new therapeutic method, recently suggested [28] [29] [30] for non-invasive 
treatment of glioma - brain cancer by a radio-frequency electric field, also opens 
new directions of understanding of glioma development. A specific task emerging 
here is whether this new medical technology is effective against invasive cells 
with a high motility, when a switching between migrating and proliferating 
phenotypes takes place. As well known, one of the main features of malignant 
brain cancer is the ability of tumor cells to invade the normal tissue away from 
the multi-cell tumor core, causing treatment failure [31] . This problem relates to 
modelling of the dynamics of cancer glial cells in heterogeneous media (as brain 
cancer is) in the presence of a radio-frequency electric field, which acts as a tumor 

1 This kind of migration-proliferation dichotomy was also found at metastatic behavior of 
breast cancer m 
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treating field (TTF) [23, 23, 30]. As reported, this transcarnial treatment by 
the low-intensity (1-3 V/cnr), intermediate-frequency (100-200kHz) alternating 
electric field, produced by electrode arrays applied to the scalp, destroys cancer 
cells that undergoing to division, while normal tissue cells are relatively not 
affectecjf|. An important result of this new technology treatment is increasing 
the survival period in twice |30pl . 

Therefore, an essential question is how the TTF affects aggressive migrating 
cells in the outer-invasive region with a low-rate of proliferation. To shed light on 
this situation, a simplified toy model for glioma treatment by the TTF has been 
suggested in 32s |33] , where a mathematical task of the migration-proliferation 
dichotomy was formulated in the framework CTRW [3U [35]. Note, that the 
simplest mathematical realization of the CTRW mechanism of the migration- 
proliferation dichotomy was introduced for a comb model [23] . In the framework 
of this toy model, it was possible to estimate the effectiveness of the TTF 
treatment in the outer-invasive region of the tumor development 132 ] . It has also 
been shown that while the TTF is highly effective in the multi-cell tumor core, 
its action is ineffective in the presence of the migration-proliferation dichotomy 
[32]. This result is mainly based on the ID consideration, where the fractal 
cancer composite is embedded in the ID space. In reality, the situation is 
much more complicated, since the fractal cancer composite in the outer-invasive 
region develops in the 3 D space. As a result of this, the TTF efficiency depends 
on the fractal dimension of the cancer composite in the outer-invasive region. 
Therefore, a more realistic model to estimate a medical effect of brain cancer 
(glioma) treatment by the RF electric field is suggested [33]. This model is based 
on a construction of a 3D comb model for the cancer cells, where the outer- 
invasive region of glioma cancer is considered as a fractal composite embedded 
in the 3D space. In the framework of this 3D model it was shown that the 
efficiency of the medical treatment by the TTF depends essentially on the mass 
fractal dimension Df r of the cancer in the outer-invasive region. 

In this paper we follow a CTRW consideration, suggested in [23] . Description 
of fractional kinetics of glioma development under the TTF treatment in the 
framework of the one dimensional (ID) and the three dimensional (3D) comb 
models show that the efficiency of the medical treatment depends essentially on 
the mass fractal dimension of the cancer in the outer invasive zone [32] [33] . The 
aim of this research is understanding both the role of the migration-proliferation 
dichotomy in fractional cancer cell transport and its influence on a therapeutic 
effect due to the TTF. 

2 An explanation of this phenomenon in the electrostatic framework is vague, since for this 
weak RF electric field of the order of lV/cm, the inter-bridge voltage between the daughter 
cells is of the order 10“ 2 3 V that is less than the voltage fluctuations related with the cell shape 
fluctuations. 

3 In the latest phase III trial study for TTF treatment for glioblastoma [§]j the TTF treat¬ 
ment has median survival of 6.6 months versus 6 months of the chemotherapy treatment. 
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2 Self-Entrapping by Fission as Fractional Mech¬ 
anism of Tumor Development 

In this section we formulate the migration-proliferation dichotomy in the frame¬ 
work of the CTRW. A simplified scheme of cell dissemination through the vessel 
network was considered by means of the following two steps [2Tij [gg] . The first 
step is a biological process of cell fission. The duration of this stage is 7/. 
The second process is cell transport itself with duration 7*. Therefore the cell 
dissemination is approximately characterized by the fission time 7/ and the 
transport time 7*. During the time scale 7/, the cells interact strongly with 
the environment and motility of the cells is vanishingly small. The duration of 
7/ could be arbitrarily large. During the second time 77, interaction between 
the cells is weak and motility of the cells leads to cell invasion, which is a very 
complex process controlled by matrix adhesion [6]. It involves several steps 
including receptor-mediated adhesion of cells to extracellular matrix (ECM), 
matrix degradation by tumor-secreted proteases (proteolysis), detachment from 
ECM adhesion sites, and active invasion into intercellular space created by pro¬ 
tease degradation. It is convenient to introduce a “jump” length X t of these 
detachments as a distance which a cell travels during the time 7*. Hence, the 
cells form an initial packet of free spreading particles, and the contribution of 
cell dissemination to the tumor development process consists of the following 
time consequences: 

Tf(l)T t (2)T f (3)... . (1) 

There are different realizations of this chain of times, due to different durations 
of Tf(i) and 7t(i), where i = 1,2,.... Therefore, one concludes that transport 
is characterized by random values T(i) which are waiting (or self-entrapping) 
times between any two successive jumps of random length X(i). This phe¬ 
nomenon is known as a continuous time random walk (CTRW) jg7]. It arises 
as a result of a sequence of independent identically distributed random waiting 
times T(i), each having the same PDF w(t), t > 0 with a mean characteristic 
time T and a sequence of independent identically distributed random jumps, 
x = X(i), each having the same PDF A(a;) with a jump length variance a 2 . It is 
worth mentioning that a cell carries its own trap, by which it is set apart from 
transport. This process of self-entrapping differs from the standard CTRW, 
where traps are external with respect to the transporting particles. The crucial 
point of the fractional transport is the power law behavior of the waiting time 
PDF 

w(t) =aT/(l + t/T) 1+a (2) 

where 0 < a < 1 and T is a characteristic time. In this case the averaged 
time is infinite. A proper explanation of eq. (0) can be the following quotation 
from Ref. [34] : “A process with the long tailed pausing time distribution would 
suffer a very sporadic behavior - long intermittencies may exist, followed by 
bursts of events. The more probable pauses between events would be short but 
occasionally very long pauses would exist. Given a long pause, there is still a 
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smaller but finite probability that an even longer one will occur. It is on this 
basis that one would not be able to measure a mean pausing time by examining 
data.” Some justification of eq. © for the fission times can be presented 
by proposing multi-time scales of self-entrapping. We can consider that self- 
entrapping for different generations of cells has different mean characteristic 
time scales, see Appendix A. One obtains that the PDF, which accounts for 
all exit events from proliferation occurring on all time scales, has the power 
law asymptotic of eq. ©■ Obtained distribution of eq. © is valid, when cell 
transport is considered on a fractional subdiffusive structure such as a comb 
model. 


3 Comb-Like Model with Proliferation 


Fractional transport of cells, namely subdiffusion, can be described in the frame¬ 
work of the comb model [38] . The comb model is an example of subdiffusive 
Id media where CTRW takes place along the x structure axis. Diffusion in 
the y direction plays the role of traps with the PDF of delay times of the form 
w(t) ~ 1/(1 + f/T) 3/2 . A special behavior of diffusion on the comb structure 
is that the displacement in the ^-direction is possible only along the structure 
axis (x-axis at y = 0). Thus, the diffusion coefficient in the ^-direction is 
D xx = Dd(y), while the diffusion coefficient in the transversal y-direction is a 
constant D yy = D 0 . A random walk on the comb structure is described by the 
distribution function P = P(x,y,t) and the current 

, N dP dP, 


The continuity equation with proliferation C(P ) yields the following Fokker- 
Planck equation 


dP 

~dt 


5{y)D 


d' 2 P 
dx 2 


Do 


d 2 P 
dy 2 


C(P), 


(3) 


where the diffusion coefficients can be related to the CTRW parameters D = 
er 2 /T. The initial condition P(X,y,0) = Po(x)5(y) is an initial distribution on 
the x axis, and the boundary conditions are taken on infinities P(t) = P'(t ) = 0 
for both the x and y coordinates. The primes denote the spatial derivatives. 

It is convenient to work with dimensionless variables and parameters. In 
the case of normal diffusion, when D x = const, the dimensionless time and 
coordinates are obtained by re-scaling with relevant combinations of the comb 
parameters D x and Dq. One obtains the following dimension variables for time 
(Do/I? 2 )t —»• t and for the coordinates D 0 x/D x —» x, Doy/D x —» y. 

We consider a possible mechanism of tumor cell proliferation. The term 
C(P) in eq. © determines the change in the total number of transporting cells 
due to proliferation at rate C. This can be considered as a linear approximation 
of a logistic population growth [39] 


C(P) = CP{ 1 - P/K ), 


(4) 
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where K is the carrying capacity of the environment (see e.g ., [40] ). It is worth 
stressing that linearization is important in the use of the powerful machinery 
of the Laplace transform. When P/K —> P < 1/2 and C = KCD\/Dq, then 
the linearization C{P) = CP is valid ©IJ. In the opposite case, when Pi >1/2 
the growth is approximated by C(P) = CP , where P = 1 — P. According to 
the migration-proliferation dichotomy in the comb model, the transporting cells 
along the x axis do not proliferate. This means that cells proliferate only if they 
have a non-zero y coordinate. Therefore, C(P) = C( 1 — d(y))P, and eq. (0 
reads in the dimensionless form 


~dt - s{y) ~d^~ a7= c(1 - <( » ))p - < 5 > 

When C > 0, eq. (0 describes cell transport with proliferation, and the PDF P 
corresponds to a low concentration of cells. In the opposite case, when C < 0, 
eq. 0 describes fractional cell transport with degradation that corresponds to 
a high cell concentration, and P exchanges for P. 

The first term in the r.li.s. of eq. © is eliminated by substitution P = e Ct F. 
Carrying out the Laplace transform F(s, x, y) = C[F{x, y, i)] and looking for the 
solution in the form F = e~^ v ^ f(x, s ), one obtains 

F{x,y,t) = C" 1 [f(x,s)exp(-i/s\y\)] . ( 6 ) 


As admitted, the true motion is in the x axis, while the y axis is an auxiliary, 
and integration over y is performed. Integrating eq. 0 with respect to the 
variable y and introducing the PDF 

/ OO 

P(x,y,t)dy, (7) 

-OO 


one obtains the following equation for F\ = e ct P\ in the Laplace space F’i(s) = 
C[Fi(t)\: 

sF 1 -d 2 J = P 0 (x)-Cf. (8) 

Integrating eq. © over y, we obtain a relation between the PDFs of the total 
number of cells Fi and transporting number of cells / in the Laplace space 

/ = F(x,y = 0, s) = (l/2) y /sF 1 (x,s). 


Substitution of this relation in eq. © yields, after the Laplace inversion, the 
Fokker-Planck equation for the distribution F. To this end, eq. © is multiplied 
by \fs and then by virtue of eq. (1C. 61) the inverse Laplace transform yields the 
following equation for F\ 

2 D^Fi - dl.Fi = -CFi , (9) 

where Dg is the fractional derivative in the Caputo form [411 ( 421 (see Appendix 
C). This equation describes fractional transport of cells with fission when C > 0 
and degradation when C < 0, where the sign of C depends on either P = e Ct F < 
1/2, or P > 1/20. 

4 Since dP = —dP, eq. J5) for P (when P < 1/2) just coincides with one for P = 1 — P 
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4 Fractional Dynamics of Untreated Cancer 


As shown, the cell fission is a source of the fractional time derivatives. This 
equation can be extended for an arbitrary fractional exponent 0 < a < 1: 
1/2 —> a. Therefore, this generalization of eq. © yields 

DgFi - adlF x = —aCF\ . (10) 


Taking into account that Dq can be expressed by the Riemann-Liouville frac¬ 
tional derivatives D^ L (see Appendix C) D g = F> c ^ L 1 D] iL and D^D 1 ^ = 1, 
we obtain another, standard, form for the fractional Fokker-Planck equation 
(FFPE) with proliferation, or degradation, 


dF 1 1 _ a d 2 Fi _ rn i- ap 

dt RL dx 2 “ c RL Fl ' 


( 11 ) 


To solve eq. (ED, we use the separation of variables [35]. We consider an 
analytical solution for the P < 1/2 using the following substitution 


Fi(x,t) = ^2,T n (t)(j> n {x). 


( 12 ) 


Therefore, a solution which corresponds to the initial condition Pq{x ), is deter¬ 
mined by the Green function G(x, t\x', 0): 

/ OO f OO /» 

dx'G(x,t\x',0)Po(x')= / dx' / dkTk(t)(f>k{x)(t> i l : (x , )Po(x'). 
-OO J— OO J 

(13) 

Here 4>k{x) is a solution of the eigenvalue problem 

d 2 4>k 


dx 2 


= \{k)4>k , 


where A (k) = k 2 is the continuous spectrum with eigenfunctions 


4>k(x) = exp 


±kx 


(14) 


The temporal eigenfunction Tk(t) is governed by the fractional equation 

T k (t) + aX c{k)D^ L a T k {t) = 0, (15) 

where A c(k) = (fc 2 +C). The solution is described by the Mittag- Leffler function 
E a (z) = E ai i(z) [13] (see Appendix C) 

T k (t)=E a [a\ c (k)t a } , (16) 

where T k ( 0) = 1, and E a (z) has the initial stretched exponent behavior 

Tk(t) ~ exp [— [a\cik)t a /T(\ + a)] (17) 

(when P > 1/2). The only difference is when P < 1/2, C > 0, while for P > 1/2 one has 
C < 0. 
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( 18 ) 


which turns over to the power law long-time asymptotics 

T k (t) ~ [T(l-a)aAc(A:)t a ] -1 . 

Using these properties of E a (z), the fractional spreading of cancer cells can be 
evaluated analytically for both initial and long-time behaviors. Substitution of 
eqs. m and in eq. (fl3l) yields the following initial time solution 


P\{x, t) oc \J nT ^ t t exp [Ct — aCt a /T(l + a)] 
x exp [—r(l + a)x 2 /4at a \ . 


(19) 


Analogously, the long-time solution is 


Pi(x, t) oc 


at a T(l — a) 


exp 


Ct — y/C\x\ 


( 20 ) 


where we take, for clarity, Pq = 5(x) for both the short and long time solutions. 
These two solutions uni) and EU1) corresponds to different scales. Solution of 
eq. (1201) describes long-time/sort-scale dynamics. When the argument in the 
exponential function is zero, it corresponds to the front of cell invasion with 
equation x ~ Iq = \[Ct. This is a so-called linear model which describes a 
solid tumor growth. In this region with x < Iq the exponential growth e Ct is 
dominant. Sundiffusion described by eq. (flUl) corresponds to the cell transport 
in the outer-invasive zone with x > Iq. When C — > 0 only this solution takes 
place. Therefore, we have the the cell spreading in the core region with x < Iq 
is due to the cell proliferation, while in the outer-invasive zone the cell motility 
is the main engine of the cell spreading. 


5 Cell Kinetics in Presence of the TTF 

Let us consider cell kinetics in the outer-invasive zone in more detail. To this 
end we consider the fractal cancer development in the presence of the TTF. This 
process can be described by fractional kinetics in the framework of the comb 
model, as well, where it is easier to draw an intelligible picture of interplay be¬ 
tween high-motility of aggressive cancer cells and the TTF in the outer-invasive 
region. Contrary to the ID comb model, in this section we extend our consider¬ 
ation of the treated cancer to the three dimensional cancer development, where 
proliferation takes place inside a fractal composite, embedded in the 3D space 
with the fractal dimension Df r < 3. 

In the 3D comb model, this anomalous diffusion can be described by the 4 D 
distribution function P = P(x,y,t), and by analogy with the ID comb model 
©, a special behavior here is the displacement in the 3D x-space at y = 0. 
The Fokker-Planck equation in the same dimensionless variables reads 

d t P = 6(y)AP + dd*P, (21) 

where d is an effective diffusion coefficient and A = X)j=i *9z ■ • 







5.1 Comb Model with Proliferation and TTF 


Obviously, cell fission/division is random in the x-space and discontinuous, 
contrary to that in the tumor core. Therefore, the outer-invasive region of the 
cancer can be reasonably considered as a random fractal set F D[r (x) = F a (x\) x 
Fp(x 2 ) x F 7 (ir 3 ) embedded in the 3 D space, as, for example, for low-grade 
astrocytomas [7], with the fractal dimension, 0 < Df r < 3 and a + /3 + 7 = Df r . 
For simplicity, we take a = j3 = 7 = Df r /3 = v. 

The effective diffusion coefficient in eq. m becomes inhomogeneous d — > 
dx( x )i where x(x) = xi x i)x( x 2)x( x 3) is a characteristic function of the fractal, 
such that xi x j) = 1 for x j G Fo fr (x) and x( x j) = 0 for Xj £ Fo fr (x), where 
Xj are the Cartesian coordinates j = 1,2,3. Now we take into account the 
influence of the TTF that affects (destroys) only quiescent cells, belonged to 
the proliferation phenotype, according to Refs. [2S1 I2S1. 3U|. Mathematically, 
this process is expressed by diffusion in the y direction with decay: 



( 22 ) 


where coefficient C defines a difference between the proliferation and the degra¬ 
dation rate. In general case, C is a random function of time and space. For 
example, it was considered as a random death rate for the random walk in 
the discrete inhomogeneous media [53] . Here we take it as a positive averaged 
constant value. 

Summarizing these arguments, mapping the glioma problem onto the AD 
comb model can be described by the following rules. (*) The dynamics of cancer 
cells takes place in the 3 D space, which is described by three x coordinates 
(x±.X 2 , £ 3 ). (H) The y axis corresponds to a supplementary coordinate that 
introduces the migration-proliferation dichotomy for the model. Therefore, at 
y = 0 the cells migrates and are not affected by the TTF. Contrarily, the cells 
with y ^ 0 proliferate and are subjected to the TTF. 

Taking this into account, one arrives at the equation of the cancer develop¬ 
ment in the presence of the TTF 



(23) 


First, we apply the Fourier transform to eq. (l23l) with respect to the Xj coordi¬ 
nates. To this end, we rewrite eq. (12311 in the form of convolution integrals. 

Therefore, as shown in |33| and in Appendix B, fractal cancer development 
in the presence of the TTF can be considered as a random fractal composite 
of cancer cells embedded in the 3D. Following coarse graining and averaging 
procedure, described in Appendix B, we arrive at the 3 D comb model that 
describes the fractal cancer development in the outer-invasive region of glioma 
in the presence of the TTF 


d t P(r,y,t) = d(y)AP(r,y,t) + [dd 2 y - C](-A) 2 f ‘ P{r,y,t ), (24) 


where P(r,y,t) is the radial function in the 3D x space r = |x|. 
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5.2 Dynamics in the Fourier-Laplace Space 

Equation (El) can be considered as a starting point of the analysis, and its 
solution will be obtained by means of the Fourier and the Laplace transforms. 
Performing the Fourier transform, constructed in the Appendix B, one obtains 
eq. El) in the Fourier space 

d t P = —k 2 5(y)P + k 3 ~ Db [ddyP - CP] . (25) 

The last term in the r.h.s. of eq. (El) is eliminated by the substitution 

P(k,y,t) = exp(-Ck 3 ~ DfI t)F(k,y,t) = e~ Ck at F(k,y,t), (26) 
where a = Df T — 3. 

The next step of the analysis is the Laplace transform in the time domain 
£[F{k,y,t )] = F(k,y,s ). 


Looking for the solution of the Laplace image in the form 

F(k, y, s) = exp[-\y\y/k a s/d]f(k,s ), (27) 


one arrives at the intermediate expression in the form of the Laplace and Fourier 
inversions 


P(r, y, t)=P k 1 


exp (—Ck a t)jC 


-l 

t 


e ~\y\y/sk<x/d 
2V sdk~ a + k 2 


(28) 


As admitted above, the y axis is the auxiliary, or supplementary coordinate, 
which determines the cell proliferating process (cell fission). Therefore to find 
the complete distribution of cancer cells in the x space, integration over y is 
performed (see Sec. W 


/ OO 

P( r ,y, t)dy. (29) 

-OO 

Both the integration over y and the inverse Laplace transform are carried out 
exactly. This, eventually, yields a solution in the form the 3 D Fourier inversion 

P{P t) = -L J” e“ ik x exp (-Ck 3 - D *t)£i ( - ± y/k'+ D *t/d) d 3 k . (30) 

Here 

. . . 1 f u a ~ 1 e u du 

£ a (-z) = — / ——— 

27 Tl J 1 u a + Z 

is the Mittag-LefHer function defined by the inverse Laplace transform with a 
corresponding deformation of the contour of the integration .431. 
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5.3 True Distributions 


Solution (l30l) is a convolution of the kernel of the TTF treatment 7 Z(z) and the 
untreated cancer distribution V(z) 

P(r,t) =1l*V . (31) 

When C = 0, which means that the TTF compensates proliferation, the solution 

1 

is described by the Mittag-Leffler function with the scaling variable z = r/t 1+D fr. 
This scaling determines the cancer cell expansion 

r ~ if (32) 


that depends essentially on the fractal dimension of the proliferation volume of 
the fractal cancer composite and reflects the migration-proliferation dichotomy 
in the outer-invasive region. Indeed, for the fractal cancer volume (or mass) 
/i(r) ~ r Dtr , the cancer development is superdiffusive when Df r < 1, while for 
Df r > 1 the latter spreads subdiffusively. This property is pure kinetic and, 
apparently, is universal for the cancer development and related to the fractal 
dimension of the cancer [?]. 

Now let us return to the convolution integral m- To avoid awkward expres¬ 
sions of integrations with the hypergeometric functions, we consider particular 
cases of the fractal dimension Df r = 2 and Df r = 1. For Df r = 2, due to the 
scaling argument, one obtains that the untreated cancer spreads subdiffusively 
(r 2 ) ~ 1 3, while for the TTF kernel we have 


H(r, t ) 


1 

\J (27r) 3 r 



e~ Ctk 0Ji(kr)dk 


1 

37r 2 (Ct) 3 


F( 3 2- 3 - 
2 7b(-,2,-, 



(33) 

Here 2-Li( a )^ c > z) is the hypergeometric function. This yields the power law 
decay of the distribution function 


K(r,t) 


(Ct)- 3 TT~ 2 
3[1 + r 2 / (Ct) 2 } 2 ' 


(34) 


This power law kernel shows that the TTF is inefficient for Df r > 1 in the 
presence of the migration proliferation dichotomy. It is tempting to calculate 
the second moment ( r 2 {t)) with the distribution P(r, t). In this case, one should 
recognize that a cutoff r = t of the Levy flights for r,l> 1 should be performed. 
This is a well known procedure [451 . used for the Levy walks. 


5.3.1 Dfr = 1 

The situation changes when Z?f r <1. In this case the TTF leads to the Brown 
exponential cutoff of the cancer spread in eq. m- For Df r = 1 the problem is 
analytically treatable. For the small argument, which corresponds (for a short 
time) to a long-scale tail of the distribution, the Mittag-Leffler function decays 
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exponentially exp ^ — Ki ^/|fc| 1+D& f j [551 S3 with the generalized transport 

coefficient K 1 = [2r(3/2)\/d] _1 . This yields the solution for the compensated 
cancer with C = 0 in the form of the hypergeometric functions like in eqs. 
and (1M1) . Following J?8], one obtains 




K iy/t. 

(2tt) 3 (1 + r 2 /Klt) 2 


Kiy/t + 
2 



(35) 


2 


This metastatic power law behavior is restricted by the Brown distribution due 
the TTF kernel 


K(r,t) = = (4icb^“ P (-icf) ' < 36 > 

The second moment is a good characteristic to show the TTF influence. One 
obtains from eq. (l35l) for the compensated cancer (r 2 ) ~ t? for r>l that 
corresponds to superdiffusion at the large scale asymptotics, and the cutoff at 
r = t is taken into account. The same calculation with the TTF kernel yields 
an effective treatment with (r 2 ) ~ t* that corresponds to the superdiffusion 
subdiffion transition due to the TTF. Obviously, that untreated cancer with 
C < 0 leads to the exponential spreading of cancer cells due to the exponential 
proliferation. 

5.4 Numerical Estimations of Eq. (1301) 

As shown in eqs. (l34l) and (1M1) analytical form of the TTF operator depends 
on fractal dimension D{ r . Since analytical estimation of eq. HU leads to awk¬ 
ward expressions of integrations with the hypergeometric functions, numerical 
procedure is performed. The results are depicted in Fig. 1 for different values 
of the fractal dimension Df r . 

As obtained, the maximal therapeutic effect takes place at D[ r = 3 that im¬ 
mediately follows from eq. HU- One should recognize that the solution for the 
compensated cancer V(r, t) is exactly the form of the interplay between the TTF 
and subdiffusion, which is the result of the migration-proliferation dichotomy. 
Another manifestation of this interplay is the fractal dimension Df r < 3 that 
leads to metastatic behavior of either the TTF kernel or compensated cancer 
solution. 


6 Conclusion 

The present study focuses on the influence of cell proliferation on transport 
properties. The mathematical formulation of this proliferation-migration di¬ 
chotomy is based on the two main stages: cell fission with the self-entrapping 
time 7/ and cell transport with durations 71- By virtue of these two time scales 
a description of tumor development is reduced to a CTRW process. A toy model 
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of cancer development is suggested by using heuristic arguments on the relation 
between tumor development and the CTRW. In this case a fractional tumor 
development becomes a well defined problem since a mathematical apparatus 
of CTRW is well established (see e.g. {331311SD H7j). The constructed model 
is a modification of a so-called comb structure [33]. An important feature of 
this consideration of cell transport in the framework of the comb model is an 
essential enhancement of anomalous transport due to proliferation. Moreover, 
we obtained that the distribution function of the fractional transport depends 
on only two parameters, namely, scaled proliferation rate C and the fractional 
exponent a, where a = 1/2 for the comb model. 

The next step is studying glioma cancer development in the outer-invasive re¬ 
gion in the presence of a tumor treating field. The model is based on a construc¬ 
tion of a 3D comb model for the cancer cell transport, where the outer-invasive 
region of glioma cancer is considered as a fractal composite, embedded in the 
3D host of the normal cell tissues. The description is performed in the four- 
dimensional (x, y) space, where the real three-dimensional x space stands for the 
description of real cancer development, while the supplementary (/-coordinate 
is introduced to described a non-Markovian process in the framework of the 
Markovian description. From the biological point of view, this corresponds to 
the migration-proliferation dichotomy of the cancer cells, where the influence of 
the TTF is considered, as well. Therefore, the kinetic equation (1331) in the (x, y) 
space is constructed by means of the coarse-graining, or embedding procedure 
inside the fractal space. This corresponds to the averaging in the 3D Fourier 
space in eq. (IB. 61) . and can be (roughly) considered as a generalization of the 
ID procedure, based on averaging extensive physical values and expressed by 
means of a smooth function over a Cantor set that, eventually, leads to fractional 
integration [351 135] . 

The efficiency of the TTF is estimated in the form of the convolution in eqs. 
m and m- This expressions describe the influence of the TTF on the cancer 
development. The efficiency of the medical treatment by the TTF depends 
essentially on the fractal dimension Df r , and the TTF is the most efficient for 
Df r = 3. But, in reality, the outer-invasive zone is a fractal composite with the 
fractal dimension Df r < 3. 

Another result relates to the spread of the compensated cancer, which is de¬ 
termined by eq. (1321) . This result reflects the migration-proliferation dichotomy, 
namely the dependence of the cancer cells spread on the fractal dimension of 
the proliferation volume. For example, the cancer development is superdiffusive 
when Df r < 1, while for Df r > 1 it spreads subdiffusively. This property is pure 
kinetic and, apparently, is universal for the cancer development with a variety 
of bio-chemical processes. Recently an experimental validation of this kind phe¬ 
nomenon has been obtained for the metastatic detaching under in vitro studying 
the breast cancer m- Important question here also is aging the treatment. Ini¬ 
tial times of the cancer development and treatment are different, and the time 
difference is unknown. Moreover, since the TTF acts on the proliferating cells 
only, the migration proliferation dichotomy leads to a particular case of a more 
general problem of aging population splitting [50]. For the present analysis this 
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general approach mm can be important for understanding the efficiency of 
the TTF. This can be an interesting issue for future studies. 

It is also worth noting that a general solution in the form of the convolution 
P(r,t ) = 7 Z-kV makes it possible to consider the compensated cancer in a more 
general framework of the fractional Fokker-Planck equation, namely [52l 

dfP = (- A)*V, 

where Caputo fractional derivative df is responsible for the migration-proliferation 
dichotomy, while q = q(Df r ) reflects the fractal dimension of the tumor devel¬ 
opment in the outer-invasive region. In our case, p = according the comb 
model construction, while q = (1 + Df r )/4 is an universal parameter, which 
determines the fractional space derivative due to the fractal dimension of the 
quiescent/proliferating cancer cells. In general case, p € (0,1), and it is de¬ 
termined by another way of the introduction of the supplementary variable y 

mug. 

In conclusion, we discuss briefly a possible direct experiment, confirming the 
existence of a fractal cancer composite in the outer-invasive region. Cancer was 
considered as a fractal composite where a random fractal inclusion of the can¬ 
cer cells F Dfr (x) = F a (x i) x Fp(x- 2 ) x F 7 (x 3 ) is embedded in the 3D space of 
normal tissue cells. Therefore, in the presence of the TTF, one can consider the 
frequency-dependent permittivities of migrating cancer cells e m and the normal 
tissue cells e n [32,. Under certain frequency of the TTF, the condition e m < e n 
can be fulfilled. These permittivities were observed in time domain dielectric 
spectroscopy in experimental studies of the static and dynamic dielectric prop¬ 
erties of normal, transformed, and malignant B- and T-lymphocytes [57] . The 
solution of the Maxwell equation for the electrostatic field in the frequency do¬ 
main yields an essential enhancement of the respond field inside the random 
fractal dielectric composite of cancer cells. Therefore, the respond electric field 
can be large enough to break the cell membrane. For example, as shown in 
137; the electric field response can be of the order of 10 4 4- 10 5 V/cm, which 
exerts the irreversible electroporation [58] due to the external TTF with am¬ 
plitude ~ lV/cm. This can be a mechanism for ablation of cancer cells, which 
effectively acts on migratory cancer cells. 

A key quantity of this cancer treatment is localization of the electroporation 
field inside the cancer. There is a straightforward analogy with nanoplasmonics 
(see e.g., [59] [60]), where the electric field enhancement is due to a so-called 
surface-plasmon resonance for a metal-dielectric composite, and localized sur¬ 
face plasmon oscillations are charge density oscillations confined to the con¬ 
ducting fractal nanostructure. The essential difference is that this biological 
cell enhancement of the electric field is not resonant, but geometrical due to 
the fractal cancer structure [613 In this connection, in vitro experiments can 
be important for further understanding the interplay between the TTF and the 
migration-proliferation dichotomy. Theoretical description of this phenomenon 
needs more realistic assumptions than those suggested here in the framework of 
the comb model. Such studies should be performed in the framework of more 
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sophisticated models of the switching between the migration and proliferation 
phenotypes 

m nans Q 2 i ng Elisa, and the next step is understanding how 
dielectric properties of cells correlate with cell motility and cell fission. Such 
experimental studies of the glioma cells can not be overestimated. 
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A Power Law PDF 

As an example, we consider that j- tli generation of self-entrapping is the Poisson 
process 

Wj(t) = t~ 1 exp (—t/Tj) 

with the characteristic time scale Tj = t j , where r = n = T is now an average 
time of cell divisions for the first generation. Therefore, following [5U 55] and 
repeating exactly the analysis of Ref. [55], we obtain, by taking into account 
events occurring on all time scales, the following distribution: 

1 — / 

w(t) = —-— ^ exp (— t/r l ) , 

3 =i 

where b < 1 is a normalization constant. Therefore, the last expression is a 
normalized sum and 

w(t/r) = Tw(t)/b — (1 — b) exp (—t/r) /b. 

Using conditions t r > 1/6, one obtains that at longer times w(t/r) = 
rw(t)/b. The last expression is equivalent to 

w(t)~ lA 1+ “, (A. 1) 


where a = ln( 6 )/ln(l/r). 

B Coarse-Graining Procedure of Fractal Cancer 
Composite 

Using the auxiliary identity 

/ Xj /* oo 

x(y)f{y)dy = -d Xj / x(y)f{y)dy 

-OO j Xj 

with the boundary conditions P(xi = ±oo) = 0 , this integration with the 
characteristic function can be carried out by means of a convolution 52, 5fi 162] . 
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Note, that 


x(y)f(y)dy= ^2 


Xj&F u ' 


f{y)S(y - Xj)dy , 


where d(y ~ %j) = y'(x) ~ \x\" 1 is a fractal density, such that on the 

finite interval (—x,x), the integral 



dy(y) ~ \x\ v 


corresponds to the fractal volume. Therefore, due to Theorem 3.1 in Ref. [51 
we have 

[ f(y)d»(y) - fq-r [ {x-y) v - 1 f{y)dy, 

Jo !» J 0 

which is defined for the finite fractal volume y-(x) = y-{xj). 

In what follows we will use the terminology and useful notations of fractional 
integration and differentiation [3511111111117]. Fractional integration of the order 
of v is defined by the operator (see Appendix C) 

-ooI x f{x) = —5-r [ f(y)(x-y) u ~ 1 dy, (B. 1) 

r(^) J oo 

1 r°° 

xC/( a; ) = p7^y f(y)(y-x) v ~ 1 dy, (B. 2) 

where 0 < v < 1 and T(^) is the Gamma function. By means of these fractional 
integration and differentiation (see Appendix C) 

w 1 _- v f{x) = d x [. 00 rj{x)\, (B. 3) 

W 1 + ~ l 'f(x) = d x [ x IU(x)\, (B- 4) 

one introduces the coarse-graining integration with the characteristic function 
in the form of the Riesz fractional derivative [56] 


x(xj)P(x,y,t) => [W 1 _~ u + W 1 + - u }P(x j ,y,t) = W 1 ~ v P(x j ,y,t ). (B. 5) 


B.0.1 Random fractal composite 

We consider a random fractal with an averaged volume, embedded in the 3D, 
which is a function of a radius only y(r) ~ r Dfr [63, 611, where r = xj . 
Therefore, the distribution function and the kernel of the fractional integration 
are the radial functions, and therefore, fractional integrations over the Cartesian 
coordinates x 3 are substituted by integrations over the radial functions. This 
averaging procedure can be performed in the Fourier space as follows 


T 


Y[W 1 - v P(r,y,t) 

3 


l[\k j \ 1 -'P({k j },y,t) => k 3 ~ D «P(k,y,t), (B. 6) 

3 
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where k = is the radius in the Fourier space and P(k, y , t) = P[P{r , y, t)] 

This averaging substitute in the 3 D Fourier space is an extension of ID embed¬ 
ding in the fractal, obtained in 62, 48], in agreement with Nigmatulin’s argu¬ 
ments on a link between fractal geometry and fractional integro-differentiation 
[48]. 49]. This is constituted in the procedure of averaging extensive physical 
values and expressed by means of a smooth function over a Cantor set that, 
eventually, leads to fractional integration [48l [49]. 

Note, that we did not use here any property of the kernel as a radial function 
that can be considered as the Riesz potential [42], as well 


1 


7 (a) 


n a '"' n |x, - x-r- ^ ii - XT-"* (y^Cx, —X'T) 

where a = — 3 and 7 (a) = 7 (Df r ) is defined by Weber’s integral 


3—a: ’ 


J zPj v {z)dz = 2 /3 T ^ 

1 of 1 
(2tt) 


u + /3 + 1 


/r 


U-p+l 


at the Fourier transform of the Riesz kernel 




,0—31 


\/k 


r a 2 Ji (rfc)dr = 7 (Df r )fc 3 Dtr . 


(B. 7) 


(B- 8 ) 


(B. 9) 


'o 


Following [42] (ch. 25), we redefine J][ -W 1 v as the fractional degree of the 
Laplace operator (—A) -Q 7 2 , namely 


Y[W 1 - v P(r,y,t) 


1 fP(r')U,dx'_ 


= )~?P(r,y,t). (B. 10) 


7 (a) J \x — a ;'| 3 

J 

This yields (see also pF5] ) 

f (-A)-t P(r, y, t) = k~ a P(k, y, t). (B. 11 ) 

Eventually, we arrive at the 3D comb model (l24l) 

d t P(r,y,t) = 5(y)AP(r,y,t) + [dd^ - C](-A)—* Jl -P(r,y,t ). (B. 12 ) 


C Fractional Integro—Differentiation 

Fractional integration of the order of a is defined by the operator 

al“f( x ) = r 7-T / f{y){x-y) a ~ l dy , (C. 1) 

B(a) J a 

where a > 0, x > a and r(z) is the Gamma function. The fractional derivative 
is the inverse operator to a /“ as a P“/( x) = a /“ a and a /“ = a D~ a . Its explicit 
form is 

aD~ a = [ X f(y)(x - y)~ l - a dy . (C. 2) 

r (-a) Ja 
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For arbitrary a > 0 this integral diverges, and as a result of a regularization 
procedure, there are two alternative definitions of a D~ a . For an integer n 
defined as n—1 < a < n, one obtains the Riemann-Liouville fractional derivative 
of the form 

aD a RL f (X) = (d n /x n )ai:~ a f(x) , (C. 3) 

and fractional derivative in the Caputo form 

a DZf(x) = a I2~ a f( n \x). (C. 4) 

When a = —oo, the resulting Weyl derivative is 

W“ = _oo D%r = _oo D a RL = . (C. 5) 


One also has -oo D^ v e x = e x This property is convenient for the Fourier trans¬ 
form 

P[W*f{x)\ = (ik) a f(k), 

where P[f(x)] = f(k). 

The Laplace transform can be obtained for Eq. (1C. 41) . If Lf(t) = f(s) is 
the Laplace transform of /(f), then 


L{DZf{t)} = s a f(s) - / (fc) (0 + )s a_1 


-1 -k 


k—0 


We also note that 


D a RL [ 1 ] = 


r 


, i>5[i] = o. 


r(l-a) 

The fractional derivative of a power function is 

„ tP- a F{/3 + l) 

RL T(/3 + 1 - a) ’ 


(C. 6) 


(C. 7) 


(C. 8) 


where /3 > — 1 and a > 0. The fractional derivative from an exponential function 
can be simply calculated as well by virtue of the Mittag-Leffler function (see 

e-9; El) : 

OO 

= E ^ • ( c - 9 ) 


k-0 


r(yfc + S) 


Therefore, from Eqs. (refB9) and (1C. 91) we have the following expression 

D a RL e xt = t a E 1 , 1 _ a (\t). (C. 10) 
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Figure 1: Dynamics of the treated cancer distribution for different values of 
the fractal dimension Df r , where from plots from 1 to 6 correspond to Df r = 
(3,2.5,2,1.5,1,0.5) 
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